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We theoretically study the optical properties of a gas of ultracold, coherently dressed three-level 
atoms in a Mott insulator phase of an optical lattice. The vacuum state, the band dispersion and 
the absorption spectrum of the polariton field can be controlled in real time by varying the am- 
plitude and the frequency of the dressing beam. In the weak dressing regime, the system shows 
unique ultra-slow light propagation properties without absorption. In the presence of a fast time 
modulation of the dressing amplitude, we predict a significant emission of photon pairs by para- 
metric amplification of the polaritonic zero-point fluctuations. Quantitative considerations on the 
experimental observability of such a dynamical Casimir effect are presented for the most promising 
atomic species and level schemes. 
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I. INTRODUCTION 

Most of the recent advances in the field of nonlinear and 
quantum optics were made possible by the development 
of novel optical media with unprecedented properties. On 
one hand, the optical response of carriers in solid-state 
materials can be controlled and enhanced by confining 
the carrier motion and/or the photon mode in suitably 
grown nanostructures P, [3, S [j, [E] ■ On the other hand, 
systems of ultracold atoms appear as very promising in 
view of all those applications which require long coher- 
ence times, e.g. quantum information processing. 

Even though the low density of an atomic gas limits the 
absolute strength of the light-matter coupling, still these 
systems have the advantage of being almost immune from 
disorder and decoupled from the environment. Further- 
more, they offer the possibility of a precise control and 
wide tunability of the system parameters in real time by 
optical and/or magnetic means. In particular, Mott [1,[3| 
(as well as band f8|]) insulator states have been real- 
ized, where a constant and integer number of atoms are 
trapped in the extremely regular potential of an optical 
lattice. Such systems constitute an almost perfect real- 
ization of the Hopfield model of resonant dielectrics . 

In the present paper we present a theoretical study of the 
classical and quantum optical properties of an atomic 
Mott insulator. Recently, the case of two-level atoms 
was investigated in [l^, [ll[ . Here we extend the Hopfield 
model to the case of a three-level system in the pres- 
ence of a coherent dressing field. The rich potential of 
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three-level configurations has already been demonstrated 
with the observation of a variety of remarkable effects, 
such as the quenching of resonant absorption by the so- 
called electromagnetically induced transparency (FIT) 
effect [H, [ll] , the light propagation at ultra-slow group 
velocities in the m/s range [H, ITsj. a.nd the coherent 
stopping and storing of light pulses Here we show 

how the peculiarities of atomic Mott insulator states can 
lead to further improvements of these experiments and, 
even more remarkably, open the way to studies of more 
subtle quantum optical effects. 

Even in its ground state, the electromagnetic (e.m) field 
possesses in fact zero-point fluctuations, whose proper- 
ties are non-trivially affected by the presence of dielec- 
tric and/or metallic bodies. One of the most celebrated 
consequence is the (static) Casimir effect, i.e. the appear- 
ance of a force between macroscopic objects due to the 
zero-point energy of the electromagnetic field [l7|, ■ In 
the last decades, this force has been the object of intense 
experimental and theoretical studies in a number of dif- 
ferent systems and its main properties can nowadays be 
considered as reasonably well understood. 

The situation is completely different for what concerns 
the so-called dynamical Casimir effect (DCE) [l9l.[20|. i.e. 
the observable radiation that is emitted by the paramet- 
ric excitation of the quantum vacuum when the bound- 
ary conditions and/or the propagation constants of the 
electromagnetic field are modulated in time on a very 
fast time scale. In spite of a wide theoretical literature 
having addressed this effect for a variety of systems and 
excitation schemes [2l|, [H, [2^ [2^ , no experimen- 

tal observation has been reported yet, mainly because of 
the difficulty of modulating the system parameters at a 
high enough speed and the presence of competing spuri- 
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ous effects. 

In the second part of this paper we show how Mott in- 
sulators of coherently dressed three-level atoms are very 
promising candidates for an experimental observation of 
the dynamical Casimir effect. As the atomic response 
to e.m. fields depends strongly on the amplitude and the 
frequency of the dressing field, a significant time modula- 
tion of the optical properties of the atomic Mott insulator 
can be induced by modulating the dressing parameters 
on a very fast time scale via standard pulse manipula- 
tion techniques 27[. On the other hand, the cleanness 
of atomic Mott insulator systems allows one to squeeze 
linewidths down to the spontaneous radiative level and 
hence to avoid those inhomogeneous broadening mecha- 
nisms that have so far limited the performances of slow- 
and stopped-light experiments. Moreover, as the dynam- 
ical Casimir radiation is collected in the optical domain 
and no relaxation processes are involved in the modula- 
tion process, no difficulties are expected to appear as a 
consequence of thermal black-body radiation or incoher- 
ent luminescence from photoexcited carriers 24 1 . 

An ab initio model is developed to confirm these expec- 
tations in a quantitative way: taking inspiration from 
recent works [2^, [2^ [s^l , we build a microscopic theory 
of the dynamical Casimir effect in atomic Mott insulators 
which explicitely includes the matter degrees of freedom. 
Thanks to the simple form of the resulting parametric 
Hamiltonian and to the relative weakness of the light- 
matter coupling constant, expressions for the emission 
intensity are obtained in closed analytical form. As ex- 
pected, the most favorable frequency region appears to 
be the middle polariton branch (the so-called dark-state 
polariton of |3l|). which shows a strong resonant coupling 
of light with the matter degrees of freedom, as well as a 
still acceptable amount of absorption losses. Advantages 
and disadvantages of the atomic Mott insulator system 
over previously studied solid state systems [H, [30| will 
then be pointed out, as well as the criteria for the choice 
of the atomic levels to be used. Quantitative estimations 
of the dynamical Casimir intensity in realistic systems 
appear as very promising in view of the experimental ob- 
servation of this still elusive effect. Generalization of the 
results to experimentally less demanding atomic states, 
e.g. Bose-condensed clouds or thermal gases is finally 
discussed. 

The structure of the paper is the following. In SecHIl 
we introduce the physical system and the model used 
for its theoretical description. In Sec lIIII we discuss in a 
systematic way the static properties of the system, such 
as the dispersion and lifetime of the elementary excita- 
tions of the system, the so-called polaritons. A calcu- 
lation of the dynamical Casimir emission in a spatially 
infinite, bulk system is presented in Sec lIVi and then 
extended to experimentally more relevant finite-size ge- 
ometries in Sec|Vl A quantitative discussion of the emis- 
sion is presented in Sec lVIl using realistic parameters of 
state-of-the-art samples. Conclusions are finally drawn 



in SeclVIB 



II. THE PHYSICAL SYSTEM AND THE 
THEORETICAL MODEL 



We consider a gas of atoms trapped in the periodic po- 
tential of a three-dimensional optical lattice with simple 
cubic geometry of lattice spacing a^. Unless otherwise 
specified, the system is assumed to be spatially homo- 
geneous with periodic boundary conditions in all three 
dimensions. The box sizes are equal to Lx^y^z, and the 
total volume of the system is V = L^LyLz- For a tight 
enough lattice potential and commensurate filling, the 
ground state of the system corresponds to an Mott insu- 
lator state, with an integer number n of atoms at each 
lattice site [1, 0] and almost negligible number fluctua- 
tions. In what follows we focus our attention on the n ^ 1 
case in which the atoms are spatially separated and do 
not interact but via the electromagnetic field. The total 
number of atoms in the system is thus N = L^LyL^ / a\ 
and the average density Uat = l/oi- The temperature 
of the system is assumed to be low enough for the zero 
temperature approximation to hold. 





(a) A scheme 



(b) Ladder scheme (c) Strongly asymmetric 
A scheme 



FIG. 1: Sketcli of the level schemes under consideration. 



The internal atomic dynamics takes place among three 
internal levels organized in either a A or a ladder struc- 
ture [12] as sketched in Fig[T](a,b) (the case of a strongly 
asymmetric A scheme of Figlljc) will be discussed in 
Sec lVIl) . The atoms are initially prepared in their in- 
ternal ground state g, which is connected to an excited 
state e by an allowed optical transition of dipole matrix 
element deg at a (bare) frequency LOe- For notational sim- 



plicity, the energy zero is set in a way to have LOg 



0. 



A coupling laser of frequency ujc dresses the atoms by 
driving the transition between the excited e level and a 
third, initially empty, state m of energy ujm- In the A 
case, the lifetime of the m state can be very long, much 
longer than the free-space radiative lifetime of the e state. 
In terms of the local amplitude £^^(1^) of the dressing 
electric field at the atomic position R, the (complex) 
Rabi frequency of the coupling is rJc'(R) = c?em Ec(^)- 
In the following we consider the case where VLc is spa- 
tially uniform |32l |: the discussion of more complex cases 
is postponed to future works. The direct transition 
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is assumed to be optically inactive dg 



0. 
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Provided the lattice potential is strong enough to fulfill 
the Lamb-Dicke condition and has the same effect on the 
atoms irrespectively of their internal state, the external 
degrees of freedom can be decoupled from the internal dy- 
namics and remain frozen in the motional ground state 
of each lattice site ISSl [SJl . 



A. The three-level Hopfield model 

A quantitative description of the many-atom system in- 
teracting with the electromagnetic field can be developed 
by generalizing the Hopfield model of a resonant dielec- 
tric [9,] to the present case of three-level atoms. In this 
approach, both the atomic electric dipole polarization 
and the radiation field are described as a collection of 
coupled harmonic oscillators. Neglecting for simplicity 
the photon polarization as well as all higher-lying pho- 
tonic bands, the vector potential operator has a simple 
"scalar" expression in terms of the photon (ph) creation 
and annihilation operators al^^ and aph.ia- 



j2'Kch ( 



ikR 



ikR 



(1) 



ke/Sz 

where the sum over k vectors is limited to the first Bril- 
louin zone (fBz) of the lattice. This approximation holds 
under the assumption of an isotropic atomic response and 
in the limit of a small lattice spacing cjea^/c <C 1 [35]. 

In our specific case of 3-level atoms, two material degrees 
of freedom are associated to every atom j, which corre- 
spond to its excitation from the g state to respectively 
the e and m states. As usual, raising and lowering op- 



erators are defined as \g)j 



)j and he J 
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for the excited e state, and analogously the a]^ ■ and 
o.m.j for the m state. In the spirit of the harmonic os- 
cillator model of @, these operators can be extended as 
creation and annihilation operators satisfying the usual 
Bose commutation rules. This bosonic description is ac- 
curate under the assumption that the probability for a 
given atom to be in an excited state is small [SOj : in this 
limit the higher-lying atomic states are not involved in 
the physics under consideration and the dynamics is tak- 
ing place mostly in the subspace spanned by the three 
|(g,e,m)) states for which the harmonic oscillator de- 
scription is exact. Generally speaking, this assumption 
is expected to be accurate as long as the number of ex- 
citations present in the system is much smaller than the 
number of atoms [3^ . 

Reabsorbing the phase of the dressing field Q.c and its 
time-dependence at u>c into the definition of the a,m,j and 
j operators, the internal dynamics of the atom is de- 
scribed by the following time-independent Hamiltonian: 

+ hVlc {al/a„^,3 + aLjOe.j) (2) 



with a real Vic and a renormalized = LOm i loc, the 
± signs referring to respectively the A and the ladder 
configuration (see Fig[T]). 

The electric-dipole coupling of the g ^ e transition to 
the transverse e.m. field [51j] 

is included by means of the standard minimal coupling re- 
placement. The Hamiltonian Q can be in fact rewritten 
in terms of the harmonic oscillator position and momen- 
tum operators as 

_ M UjI ^2 , 1 p2 , 

Mujf^ ^2 , J_ p2 , 

2 "'J 2A/ "'^ 



PejPm.j (3) 



where 



Pe,j — i 



(4) 



(5) 



Analogous expressions hold for Xm.j and Pmj, where 
dej is replaced by dmj and tOe by (D^. The minimal 
coupling Hamiltonian is then obtained by replacing P^j 

with Pej — qA{'Rj)/c in while Pm.j is not affected 
because the g m transition is optically forbidden. Rj 
is here the position of the j atom. 

Note that the "charge" q and "mass" M parameters that 
appear in the harmonic oscillator model do not have a 
physical meaning per se, but are model parameters that 
are to be chosen in a way to reproduce the physics of 
the system under examination. By comparing the dipole 
moment of the atomic transition to the matrix element 
between the ground and the first excited state of the 
harmonic oscillator model, one finds that they have to 
satisfy the relation deg — y/hq"^ /2MuJe- All observable 
physical quantities will in fact involve only deg and not 
the q and M separately. 

To take full advantage of the translational symmetry of 
the system, it is useful to introduce the collective atomic 
operators 
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a 



(e,m) ,k 



Ed] X • e 



ikRj 



(6) 



which create a delocalized atomic excitation with a 
wavevector k belonging to the first Brillouin zone of 
the lattice. Analogously to their localized counterparts 
a(e,m),j and a]^e,Tn),j' a(e,m),k and dj^ satisfy Bose 
commutation rules. 

Straightforward manipulations lead to the final form of 
the light-matter Hamiltonian: 



H 



H, 



ai.k 



ini ,kj 



(7) 
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where 

HatM 



■ k + flp/i -k) (op'' -k 

(am,k - aL,-k) (V 



1 



(8) 
(9) 



+ hDk [ciphM + V,-k j (ap'':-k + V,k 



2^ 



k()o) 



All terms consist of quadratic forms in the creation and 
annihilation operators of the electromagnetic or the mat- 
ter polarization fields. The first term -ffp^^k is the free 
e.m. field Hamiltonian. The second term ffat.k describes 
the internal dynamics of the dressed atoms. The three 
lines of the third term -ffi„t_k respectively account for (i) 
the dipole coupling of the g e transition to the e.m. 
field, (ii) the photon renormalization due to the squared 
vector potential term, and (iii) the coupling between the 
photon quantum field and the m excitation as a result of 
the dressing field 52]. 

The coupling constant Ck is equal to 



hck 



den , 



(11) 



and Dj, = C^/wg. As expected, these expressions involve 
the model parameters q and M only via their physical 
combination deg. 

Introducing the vector dk of bosonic operators: 

ttk = (aph,k, ae,k, am,k, ap;i,_k' '^e.-k' "m,-k) ' (l^) 

and the Bogoliubov metric 77 ~ diag[l, 1,1,— 1,-1,— 1], 
the Hamiltonian ([7|) can be recast in a simple matricial 
form: 



h 



(13) 



in terms of a 6 x 6 Hamiltonian matrix Tik of the form: 

{-%-% )■ 

where ICk and Qk are 3x3 matrices. The constant Eq 
fixes the energy zero: as it has no consequences in what 
follows, it will be neglected from now on. 

The Hermitian matrix 



ck + 2Dk iCk iflcCk/uJe 

A^k = I —iCk UJe 



(15) 



takes into account the free field, the internal atomic dy- 
namics including the dressing beam, as well as the light- 
matter interaction terms at the level of the so-called Ro- 
tating Wave Approximation (RWA): whenever a radia- 
tive photon is absorbed (emitted), an atomic excitation 
is created (destroyed) at its place [s^]- 
The symmetric matrix 



2Dk -iCk -iVLcCk/i^e 
Qk = I -iCk 

-iVlcCkluJe 



(16) 



corresponds instead to those additional terms which de- 
scribe anti-RWA, off-shell processes where a photon and 
an atomic excitation are simultaneously destroyed or cre- 
ated. 

The relative importance of the RWA /Ck and the anti- 
RWA Qk terms is quantified by the ratio C/uje of the 
radiation-matter coupling strength C = Ck^^^/c and the 
excitation frequency oje- For most atomic systems of ac- 
tual experimental interest, this parameter is generally 
quite small: as a simplest example, consider the D2 line 
of ^^Rb atoms at Ag = 2tt c/uj^ ~ 780 nm. The electric 
dipole moment of the transition is d^g — 'i.2eaBohr [39j 
and a typical value of lattice spacing is — 300 nm. 
For a unit filling factor n — 1, the light matter cou- 
pling parameter is then C/uJe — 1-7 • 10~*. Although the 
condition C/uje ^ 1 rules out the possibility of observing 
the so-called ultra-strong coupling regime [2^, [13] in such 
dilute atomic systems, the anti-RWA terms in the Hamil- 
tonian can still have interesting observable consequences 
as we shall see in what follows. 



III. STATIONARY STATE: GROUND STATE 
AND POLARITON EXCITATIONS 



We begin the study of the optical properties of the system 
from the simplest case where the dressing parameters ujc 
and Qc are kept fixed in time. For each value of them, 
the quadratic structure of the Hopficld Hamiltonian (fT5|) 
guarantees that this can be set into the canonical form: 



H 



(17) 



by means of a Hopfield-Bogoliubov transformation 
As in (fT3|) . the constant Eq is the zero-point energy and 
will be neglected from now on. For each wavevector k, 
the frequencies uJr,k of the elementary excitations are 
given by the eigenvalues corresponding to the positive- 
norm eigenvectors of the Hopfield-Bogoliubov matrix 



Mk = {-HkY 



(18) 



In the system under consideration here, the elementary 
excitations are the lower (r = LP), middle (r = MP, 
often also called dark-state polariton, e.g. in ^3]) and 
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upper (r = UP) polariton modes. All these modes are 
linear superposition of light and matter excitations. The 
polaritonic annihilation operators pr^k can be written in 
terms of the eigenvectors Wr^k as: 



E 

A=l 



(19) 



The index A runs over the six components of the eigen- 
vector Wr,k of the Hopfield-Bogoliubov matrix p8|) and 
of the operator vector defined in p^ . An analogous 
expression holds for the creation operators p\^. 
Grouping the p^.k's in the operator vector 

TTk = (pLRk,PA^p,k,P[/p,k,p[p_^,pj^,^p_^^p|^p_^)^, 

(20) 

the transformation ([TO]) to the polaritonic basis can be 
cast in the simple matricial form TTk = WkAk; for in- 
stance, the column vectors obtained as the transposed 
of the first three lines of Wk correspond to the iw^.k eigen- 
vectors (HSl), W{i,2,3},k = W{Lp^MP,UP}.Vi- 

The orthonormality condition of the eigenvectors in the 
Bogoliubov metric rj corresponds to the 77-unitary condi- 
tion 



and guarantees that the operators p^.k satisfy the stan- 
dard Bose commutation rules. In the basis, the Hamil- 
tonian matrix has the simple diagonal form: 

^k = WkTik = diag[wLP,fc,t^MP,fe,t^;7P,fc, 

— '^LP,k,—'^MP,k,—^UP,k]- (22) 

In view of the following developments, it is useful to give 
the explicit form of the r = {LP, MP, C/P} polariton 
operators in terms of the photonic {ph) and matter (e, m) 
excitation ones: 

Pr,k = Mr,k °-phM + ^r.k ^e.k + w"k Om.k 

+ <'k «^.,-k + <k «l^k + <k «L,-k (23) 
as well as the inverse transformation (j = {ph, e, m}): 



aj,k = ^LP.kPLPM 



' "^MP.kPMPM 



^UPkPUPM 



^ip,kPip,-k ^MP,kP\lP,-k ^UP,kPuP,-k- (^4) 

The u and v Hopfield coefhcients characterize respec- 
tively the normal and anomalous weights of the different 
ph, e, m components of the LP, M P,UP polaritons. 



A. The polariton vacuum 

The vacuum state of the system corresponds to the 
ground state |G) of the Hamiltonian (fTT)) . and is defined 
by the vacuum condition 



PrM |G) = 



(25) 



for all polariton modes r = {LP, MP, UP}. 

As both annihilation and creation operators are in- 
volved in the Bogoliubov transformation (|24p , the ground 
state |G) corresponds in the original a(p^ e,m) basis to a 
squeezed vacuum state with a non-vanishing expectation 
value of the photon and atomic excitation numbers: 



= E 

r={LP,MP,UP} 



(ph,e,m) i2 



(26) 



As the atomic systems under consideration here are far 
from the ultra-strong coupling regime [1,[2^,[33], an accu- 
rate estimation of N9, > can be obtained by means of 
perturbation theory in the light-matter coupling strength 
C /uje ^ 1. The dressing amplitude fJc is assumed to be 
at most of the order of C. 

The zeroth order approximation j, of the eigenvector 
can be obtained by diagonalizing the block diagonal ze- 
roth order Hopfield-Bogoliubov matrix 



Ml 





-/Ck 



(27) 



(21) This provides the zeroth order eigenvalues and 



eigenvectors W;^ k- These latter have the form 

<k = K!k°.<l""k°. 0,0,0)^ (A = 1,2, 3), (28) 

WA.k = (0, 0, 0, <_3.k, UA-s.k) = 4, 5, 6). 

(29) 

Up to this level of approximation the virtual occupation 
(|26p is then rigorously vanishing. 

When evaluating the first order correction, attention has 
to be paid to the non-positive nature of the 7y-metric: 



'A.k 



E 

A' = {1...6} 



w°/. '7'57WfeW^,k 



'A'.ki 



(30) 



■'A.k 



■'A'.k 



where the sign ex is +1 for A = 1,2,3 and —1 for 
A = 4, 5, 6. The perturbation matrix 6Mk = Mk — 
is of order C /uJe and has non-zero entries only in the off- 
diagonal 3x3 blocks. Keeping in (|5D|) only the lowest 
order terms in C/uJe, one gets to the first order correc- 
tions: 



E 



-iCk 



r.k 



e,0 ph,0\ j,0* 



(31) 

where the sum runs over the three polariton branches 
r' = {LP,MP,UP}. Note that the off-diagonal terms 
due to the squared vector potential give a contribution 
to ([3T|) whose amplitude is of the order of Dk/cOe — 
[Ck/iOeY and have been therefore neglected. The same 
for the off-diagonal terms due to the dressing field, whose 
contribution is of the order of Ck^c/^l- 
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which leads to anomalous amplitudes 



FIG. 2: Hopfield it and v coefficients (a,b), the group velocity 
(c) and the absorption coefficient (d) of the middle polariton 
(MP) as a function of dressing amplitude flc in the resonant 
{uje = Cjm) case. On the scale of the figure, the analytical 
approximations (|40|l . (|35II37[) and p9|l are undistinguishable 
from the exact calculations. Vertical lines indicate the dress- 
ing amplitude values used in the next figures. The system 
parameters correspond to the D2 line of a Mott insulator sys- 
tem of ^'^Rb atoms with filling factor n = 1 in a, ul ~ 300 nm 
lattice: C/ue = 1.7 • 10"", Xe = 2-K/ke = = 780 nm. 



The virtual population in the ground state is then of 
the order of {Ck/i^eY'- exception made for a small re- 
gion around k = where Ck diverges, this population 
is therefore very small for all polariton branches. This 
fact provides an a posteriori justification of the use of 
perturbation theory. 

In the most significant resonance region k kg = uJe/c, 
the frequency denominator of (|3ip can be approximated 
by 20^6 in a sort of degenerate polariton approximation: 



ph,l 



e,l 



iC 
iC 

2uJe 



phfl 
'■r,k ' 



(32) 

(33) 
(34) 



Compact formulas for the fully resonant point Ue = Corm 
k = ke oi the MP are immediately obtained by inserting 
in (|32j|34p the explicit form of the RWA eigenvectors of 
the zeroth order Hopfield-Bogoliubov matrix l[27|) : 



ph,0 
''^MP,k^ 



e,0 
^MP,k^ 

m,0 
'^MP,k, 



((72+^(2^)1/2^ 



iC 

"(C2+r!2,)l/2 



(35) 
(36) 
(37) 



V 



ph,l 
'^MP,kc 

e,l 

MP,kc 



^MP.kc 



iCQc 



2We(C2 



f72^)V2- 



(38) 

(39) 



The accurateness of these analytical approximations is 
visible in Figl2l^a,b) where the Hopfield u and v coeffi- 
cients are plotted for the fully resonant k = k^ point 
of the MP as a function of the dressing amplitude: on 
the scale of the figure the analytical approximations are 
undistinguishable from the exact calculations. 

As long as the system parameters are kept constant in 
time, the virtual populations (|26p are intrinsecally bound 
to the system ground state and can not be revealed 
by a standard photodetector based on absorption pro- 
cesses [2^, [s^]. On the other hand, the dependence of 
the anomalous amplitude ((39|) on flc [see Figl2jb)] sug- 
gests that the zero-point fluctuations in the ground state 
can be externally controlled by varying Qc in time. In 
particular, if the time-modulation of is sufficiently 
fast, the system is not able to adiabatically follow the 
time-dependent vacuum state. As a consequence, real 
polaritons are created in the system and then emitted as 
radiative photons into the surrounding free space where 
they can be detected by a photodetector. This will be 
the subject of the next Secs lIVlTVll 



B. Polariton dispersion 

Examples of polariton dispersion are shown in Fig[3] for 
the most significant cases. The shape of the three po- 
lariton branches changes in a substantial way depending 
on the dressing parameters ujc and fip, which provides 
a simple way to externally vary the optical properties 
of the system in real time. An application of atomic 
Mott insulators as dynamic photonic structures was in- 
deed proposed in [3a |. 

As long as linear optical properties are considered, it 
is important to note that the predictions of the quan- 
tum model are indistinguishable from the solution of 
the Maxwell equations including the semiclassical expres- 
sion for the local dielectric polarizability of three-level 
atoms In particular, this is the case of the band 

dispersions shown in FigO 

As LOe and ujm have no spatial dispersion, the three bands 
have no spectral overlap and are separated by energy 
gaps in which the radiation can not propagate. For the 
systems under consideration here, the gaps between the 
bands are however very narrow, on the order of /uje, 
and therefore almost invisible on the scale of the figure. 

When the dressing is on resonance with the m — *■ e tran- 
sition (we = (I'm), the atomic resonance is split into a 
symmetric Autler-Townes doublet at tOe ± and the 
photonic mode anticrosses each of the two components. 



7 



4x10 



S 



-4x10 




4x10^ 



0.002 



-0.001 



0.001 



4x10 



S 
S 



-4x10 



-0.002 -0.001 



0.001 



^ 8x10 — 



3 
S 



4x10 



-4x10 



0.002 -0.001 0.001 

k/k - 1 



0.002 



-(C) 




1 




/ ^ — 












MPy/ 




1 


1 


1 



0.002 




0.002 



— 



-4x10 



5 Lu 




T 1 1:1 I I I I I 



I I I I miiiil I I I I I I 



lO '" lo ' 10 10 ' 10 10 ' 10 * 10"^ 10"^ 10 ' 10° 



4x10 



-4x10 



I I Kllllj 

— 




i 'l'mi I I I I I I I I 



lo " 10 ' 10 ' 10 ' lo ' 10 ' 10 ^* 10 ' 10 ' 10 ' lo" 



iiiiiiiq iiiiiiii iiiiiiiii iiiiiiiii iiiiiiiii iiiiiiiii 




I I I I I I 



lO '" 10 ' 10 10 ' 10 10 ' lO"* 10"^ 10"^ 10 ' 
abs , , 

K /k 



10 



FIG. 3: Dispersion relation (a,c,e) and absorption coefficient (b,d,f) of the three polariton modes. Same system parameters as 
in Fig[2] Blaclc, red and green lines refer to the lower polariton (LP), the middle polariton (MP), the upper polariton (UP), 
respectively. The blue dashed line in (a,c,e) is the free photon dispersion lo = ck. Panels (a,b): resonant uie ~ t^m, weak dressing 
Qc/C = 1.86 • 10"'' <C 1 case. Insets in (a,b): magnified views of the MP. Panels (c,d): resonant uje ~ Cjm, strong dressing 
Qc/C = 2 case. Panels (e,f): non-resonant Cjrn — cue = 5 C , strong dressing flc/C = 2 case. In the absorption (b,d,f) panels: 
A level scheme with 7e/27r = 6 MHz and 7m/27r — 10 Hz (solid lines); ladder level scheme with exchanged 7e/27r — 10 Hz and 
7„i/27r — 6 MHz values (dotted lines). 



Two sub-cases are to be distinguished: for a strong dress- 
ing flc ^ C [panel (c)], the two anticrossings are al- 
most completely separated and have a half-width equal 
to in between the two anticrossings (i.e. in the 

neighborhood of We), the middle polariton (MP) almost 
coincides with the photon branch and has a group veloc- 
ity v^p f. = dujM p,k/dk close to c. On the other hand, for 
a weak dressing flc C, the two anticrossings overlap, 
which results in a strong mutual distortion [panel (a)]: 
the MP dispersion is strongly flattened, and its group 
velocity becomes orders of magnitude slower than c [in- 
set of panel (a)]. 

An approximate, yet quantitatively very accurate expres- 
sion for the group velocity v ^ of the MP around res- 
onance k = fee = tJe/c is easily obtained from the MP 
eigenvector (I35II37P of the RWA Hopfield-Bogoliubov ma- 
trix (123: 



gr 



_ I ph,0 
— C \Uj^pj^^ 



(40) 



As one can see in Figl^Ja), arbitrarily slow values of 
^MP k '"^'^ obtained by simply reducing the ampli- 
tude Vic of the dressing field: thanks to the trapping 
of atoms at lattice sites, no lower bound to the group 
velocity appears as a consequence of the atomic recoil 
in absorbing and emitting photons [l^l- For the specific 
case of Rb atoms considered in the figures, a quite con- 
servative value ric/27r — 12 MHz (i.e. twice the radiative 
lifetime of the e state) already leads to v^J =11 m/s. 

In the case of an non resonant dressing [panel (e)] , the os- 
cillator strength of the optical transition is shared in an 
asymmetric way by the two components of the Autler- 
Townes doublet. This fact is responsible for the anti- 
crossing to be wider for the component having a larger e 
state weight (the lower one in the figure) . Because of the 
optical Stark effect, the position of this main anticrossing 
is slightly shifted by /S.ujose = ^'c/i'^e — ^m) from its 
bare position at LOe- 
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C. Effect of losses 

Losses from both the e and the m states are responsible 
for the finite hfetime of polaritonic excitations. For typ- 
ical systems, both decay rates 7(e,m) &re much smaller 
than the radiation-matter coupling C and the energy 
splitting ojq^k — ^q' ,k between different polariton bands 
q ^ q' a,t 'A given wave vector k. In optics, this regime 
goes often under the name of strong coupling regime [2, 13| 
and is characterized by losses not being able to effectively 
mix the different polariton branches, which retain their 
individuality. 

The Fermi golden rule then provides an accurate predic- 
tion for the decay rate 7r,k of the plane wave polariton 
state of wave vector k on the r branch: 



m |2. 



(41) 



physically, this decay rate is the relevant one in a light 
stopping experiment where light is stored in the system 
during a macroscopically long time [T6l,l3l|. In a propaga- 
tion configuration (i.e. at given frequency), the relevant 
quantity is instead the absorption length [4l| : 



V,k 



V,k 



hr 



(42) 



The relative value of the 7(e.m) loss rates depends on the 
specific level scheme under consideration. 

In a A configuration, 7^ is mostly given by non-radiative 
effects, and generally has a very small value. As no in- 
trinsic effect is expected to significantly contribute to 7™ 
in Mott insulator states, it can in principle be suppressed 
to arbitrarily small values by means of a careful experi- 
mental setup. On the other hand, radiative decay from 
the e state into the m state can take place by sponta- 
neous emission of a photon into a mode different from 
the coherent dressing one. The rate 7e of this process is 
equal to the one of an isolated atom in free space [ssl ]. 
and can be reduced only by choosing an optically weak 
TO ^ e transition to dress the system. 

In a ladder configuration, the roles of 7e and 7™ are ex- 
changed. Spontaneous emission processes can now occur 
on the TO — !■ e transition giving rise to a significant 7^ 
value. On the other hand, provided the g e transition 
is closed and no other decay channel is available for the e 
state, no spontaneous process can contribute to 7e, which 
can then be made arbitrarily small by means of a careful 
design of the experimental set-up. 

This remarkable fact was discussed at length in the pio- 
neering paper [9^, where it was pointed out that polari- 
tons in a lattice of two-level atoms are not subject to 
dissipation and can propagate with a slow group velocity 
along macroscopic distances. This fact originates from 
the rigidity of the lattice that freezes the atomic motion 
in the ground state of each site and prevent atoms from 
incoherently scattering the light. As recently suggested 
in [10j,|Tl|, atomic Mott insulators then constitue an ideal 
realization of this two-level Hopfield model. 



Plots of the prediction (|42|) for the absorption coef- 
ficient for the three polariton branches are shown in 
Figsl3^b,d,f) for respectively the A (solid lines) and the 
ladder (dashed lines) three-level schemes [sS]. Again, 
these curves are in perfect agreement with the solution of 
the Maxwell equations using the semiclassical expression 
for the dielectric polarizability of three- level atoms fl^. 
As expected, absorption is strongly peaked in the anti- 
crossing regions where polaritonic bands have the largest 
weight of matter excitations and the slowest group veloc- 
ity. This, in both cases of A and ladder configuration. 

The only exception is the dip that is visible exactly on 
resonance with the two-photon Raman transition lu = ujm 
in the case of a A configuration. In optics, this ef- 
fect goes under the name of Electromagnetic Induced 
Transparency (KIT) effect [ll[l3: in the vicinity of the 
resonance, quantum interference suppresses the weight 
of e excitation in the MP mode. Consequently, the 
absorption rate is quenched to the non-radiative one 
iMP.k^ — Im ^ 7e- This cffect is most dramatic in the 
case of a weak dressing fie ^ C [see in particular the 
left inset of FigjJJb)] , where the minimum of the spatial 
absorption coefficient K^fPfc around resonance remains 
remarkably deep in spite of the very slow group velocity 

This allows for MP polaritons to propagate for macro- 
scopic distances without being appreciably absorbed [55| : 
as one can see in Figl^l^b), a group velocity as slow 
as v^jpj, = 15m/s still corresponds to an absorption 
length as large as ^MP.fce — 16 cm. 

To conclude the section, it is useful to summarize the 
main advantages of using a Mott insulator for slow-light 
experiments: 

1. As shown in the pioneering paper Q, slow light 
propagation can be observed in a lattice of two- 
level atoms almost without absorption. 

2. Systems of three-level atoms allow for a real-time 
external control of the optical properties (polariton 
dispersion, group velocity, zero-point fiuctuations) 
by simply varying the frequency ujc and the ampli- 
tude flc of the dressing beam. 

3. As the atoms interact with each other only via the 
electromagnetic field, no intrinsic decoherence ef- 
fect can contribute to the 7m rate in a A configu- 
ration nor to the 7e one in a ladder configuration. 

4. The presence of a single atom at each lattice site 
eliminates the inhomogeneous broadening of the 
transitions that originates e.g. from the spatially 
varying density profile of a Bose-Einstein conden- 
sate. 

5. The trapping of atoms at the lattice sites eliminates 
the lower bound to the group velocity that atomic 
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recoil would impose in the case of a homogeneous, 
untrapped gas 



IV. DYNAMICAL CASIMIR EMISSION IN THE 
PRESENCE OF A TIME MODULATION 



When the dressing parameters are modulated in time, 
the quadratic form of the Hamiltonian is preserved, yet 
with a time-depedent Hamiltonian matrix Hk(i) = + 
SHk{t). At each time t, a Bogoliubov transformation 
diagonalizing the istantaneous Hamiltonian can be still 
found, but the transformation matrix Wk(i), as well as 
the polariton bands iOr^k and the expression of the polari- 
ton operators pr,k in terms of the original aj,k ones are 
now varying in time, as well as the vacuum state \G{t)) 
of the system. While for slow modulations the system is 
able to adiabatically follow the istantaneous ground state 
\G{t)), excitations are created in the case of a faster mod- 
ulation. The study of the properties and the intensity 
of this dynamical Casimir emission is the subject of the 
present and the next sections. 

Our strategy is to reduce the problem to a simple and 
tractable parametric Hamiltonian to which the standard 
tools of quantum optics can be applied. In particular, we 
shall concentrate our attention on the simplest case of 
a weak time-modulation ||(5Hk|l ^ ll'^^kll for which per- 
turbation theory can be used to obtain analytical pre- 
dictions: most among the perturbations that one can 
envisage to apply to the atomic system largely fulfill in 
fact this condition. A discussion of the physics beyond 
perturbation theory can be found in [2l|, [2^, IH, H^l 
Let Wk be the Bogoliubov transformation diagonalizing 
the unperturbed Hamiltonian matrix Tik- In general, the 
perturbation Hamiltonian 



(43) 



is not diagonal in this basis. In terms of polariton cre- 
ation and annihilation operators, STi.'-^ introduces terms 
of two kinds [4^. The diagonal 3x3 blocks corresponds 
to terms of the form pj^j^Ps.k, which are responsible for 
a renormalization of the polariton energies (for r =^ s) 
and for the occurrence of interbranch transitions (for 
r s) which transfer already existing polaritons from 
one branch to another. At the lowest order, these terms 
have no effect on the polaritonic vacuum state and will 
therefore not be considered in what follows. 

The off-diagonal 3x3 blocks of 7i[. are more inter- 
esting in the present context, as they correspond to 

terms of the forms pj^k^'l-^k ^^'^ Pr,iiPs,-ii, which re- 
spectively create or destroy pairs of polaritons in the 
opposite ±k momentum states of the r and s branches 
(r, s = {LP, MP, UP}). In particular, they account for 
the creation of correlated pairs of polaritons out of the 
vacuum state of the system via parametric amplification 
of the zero-point quantum fluctuations. 



This emission of radiation is an example of the still unob- 
served dynamical Casimir effect (DCE) ^9|]. The most 
celebrated example of DCE is predicted for a metallic 
cavity whose length is varied in time by means of a me- 
chanical motion of its mirrors 12 1| . Another possibility 



consists of varying the effective len gth of the cavity by 
modulating the mirror conductivity [2^, [2^, by mimick- 
ing moving mirrors via a suitably chosen x^^^ non-linear 
optical element j25j . or by varing the bulk refractive in- 
dex of the cavity material [12, E^] . 

A generalization of this latter scheme is the subject of the 
present paper: the modulation of the dielectric properties 
of the atomic medium is created by varying in time the 
dressing field amplitude. Using the formalism developed 
in the previous sections, we will be able to go beyond the 
non-dispersive dielectric approximation made by most of 
the existing works, and fully include the resonant dynam- 
ics of matter excitations. Using this microscopic model, 
accurate predictions will be obtained also for the fre- 
quency window in the neighborhood of the atomic res- 
onance where light-matter coupling is the strongest as 
well as the expected intensity of the dynamical Casimir 



emission. 



A. Polariton emission rate in a bulk geometry 

Depending on the frequency spectrum of the time- 
modulation 6Ti.k{t), polaritons can be emitted in any 
momentum state: a monochromatic oscillation at a fre- 



quency u! of the form SHuit) — SHi^ (e" 



0, 



fact able to resonantly create pairs of polaritons in the 
r, s branches at wave vectors ±k fulfilling the parametric 
resonance condition 



Wr,k 



-k = ^■ 



(44) 



Note that an analogous condition was recently obtained 
for the parametric emission of phonons in trapped atomic 
Bose-Einstein condensate in an optical lattice when the 
lattice potential is modulated in time |i45i] . 

Starting from the polaritonic vacuum as initial state, and 
limiting ourselves to lowest-order effects in the modula- 
tion amplitude, we can safely neglect the pl^ ]^PsM terms 
that arise from the modulation and rewrite the time- 
dependent system Hamiltonian (|13p in the standard para- 
metric form 14611: 



H = H 



H 



SH 

h 
+ 2 



El 



K-s.kC 



h.c. 



(45) 



From this Hamiltonian the polariton creation rate is de- 
termined by means of the Fermi golden rule. As final 
states, one has to consider pairs of polaritons created in 
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the r, s = {LP, MP, UP} branches at wavevectors ±k: 



dNr 



dt 



k 



(46) 

In the weak modulation regime (i.e. far below any para- 
metric oscillation threshold), this approach is equivalent 
to other ones based e.g. on the input-output formalism; 
including all neglected terms is instead crucial if one is 
interested in the peculiar pulse shaping and frequency 
up-conversion effects discussed at length in [2l], [30| • 

Replacing as usual the the sum over k vectors with an 
integral, the total creation rate per unit volume into the 
r, s branches reads: 



dNrs 

dtdV 



27r(v3^.+v%) 



I ^ rs,fe I 1 



(47) 



where the r (s) polariton is assumed to be emitted in the 
fcz > [kz < 0) half-space, k is the wavevector value 
at which the resonant condition is satisfied for the 



r, s branches under examination, v^^ ^ are the group 

velocities, and the matrix element Vrs.k of the process is 
given by: 



Vrs,k = Wk<5HkWk' 



r,s+3 ' 



(48) 



here, the q,s = {LP, AIP,UP} branches correspond to 
respectively {1,2,3}. Perturbative expressions for the 
matrix element P5|) to the leading order in C/wg will be 
given in the next subsection Sec lIVBI 

It is worth noticing that a finite emission intensity 
{pI i^Pr,k) 7^ is obtained while the classical amplitudes 
remain strictly zero PrM = (Pr,k) = during the time- 
modulation. The dynamical Casimir emission is in fact 
a purely quantum effect due to the parametric amplifica- 
tion of the zero-point fluctuations of the polariton held. 

Eas. ([T7)l and are the central result of the present 
section; they quantify the polariton emission in an ideal- 
ized, spatially infinite system and will represent the cen- 
tral building block in the study of experimentally relevant 
finite-size geometries that we shall perform in Sec|Vl A 
quantitative discussion of the actual value of the emission 
rate for realistic values of the system parameters will be 
given in SeclVIl 



B. Approximate analytical expression of the 
matrix element 



A simple, yet accurate estimation of the matrix ele- 
ment (j48p can be obtained by means of the perturba- 
tive approximation of the eigenvectors Wr^k discussed in 
sec lIII Al As a most significant example, we consider a 
periodic modulation of the coupling amplitude Qc of the 
form: 



nc{t) = nc + snc {e"^' + e'"*) 



(49) 



while its frequency is kept constant [56| . Modulations 
of this kind were at the heart of recent light-stopping 
experiments [l^, albeit on a much slower time-scale: in 
that case, the modulation had in fact to be adiabatic 
enough not to induce inter-branch transitions. 

Inserting the perturbative result (|32ll34p into the Wk 
transformation matrix, the following expression for the 
matrix element is found, which is valid throughout the 
whole resonance region: 



iC 



■ ^s,k^r,k 



(50) 



As expected, Vrs,k is proportional to the amplitude of the 
dressing modulation and to C/2uJe, i.e. the ratio between 
the amplitude of the antiresonant light-matter coupling 
and the energy associated to the creation of a pair of 
excitations. The presence in (|50|) of both photonic and 
m excitation Hopfield coefficients means that a mixing 
between light and matter modes is necessary to obtain a 
sizeable emission. A most favourable region is therefore 
the fully resonant point where LUe — LUm and the modu- 
lation is driven at a frequency uj — 2uje. 
An analytical estimation of Vrs,k^ with r = s = MP is 
readily obtained for this case by inserting in ([50)1 the an- 
alytical eigenvector p5ll37p of the zeroth order Hopfield- 
Bogoliubov matrix (|27|) : 



Vj 



MP.kc 



(J2 



SQc- 



(51) 



In the slow light regime flc <C C, the matrix ele- 
ment VMP,k^ grows as flcS^c/'^e, while it goes as 

S^c/i^JuJe) for > C. 
Using the expression (HOI) for the group velocity, it is 
immediate to see that for a given value of the relative 
modulation amplitude 5flc/^Ci the emission intensity 
(|T7)) starts proportionally to fl^ for small values ^ C 
and then saturates to a finite value for Oc ^ C. We will 
come back on these issues in the quantitative discussion 
of SecED 



The main result ([50]) of the present section fully includes 
the dispersion of polaritons, which is crucial to correctly 
describe the region around resonance. This, in contrast 
with previous works [43j where the Casimir emission was 
studied for a dispersionless dielectric medium with a pe- 
riodically oscillating dielectric constant 



6(t) 



6e (e^' 



(52) 



Also in this Hamiltonian of the parametric form 

(|45p can however be written, with the matrix element 



Vk^^.Se. 
4e 



(53) 



As a single photonic branch is considered, no r, s indices 
are needed. A quantitative comparison of the two cases 
will be made in the final section Sec lVIl 
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EMISSION RATE FROM FINITE-SIZE 

SYSTEMS 



In order to obtain a quantitative prediction for an actual 
experimental set-up, one has to go beyond the idealized 
infinite geometry system considered so far, and study the 
more realistic case where dynamical Casimir light is gen- 
erated in a spatially finite system and then revealed by 
a detector located in the external vacuum. 

The first step is to rewrite the parametric Hamiltonian 
P5|l for a bulk system in a local form. Expressing the 
polariton operators PrM in terms of the real-space ones 



1 



p^,k = I due *''^Vr(x), 



(54) 



the parametric Hamiltonian can be rewritten as: 
SH^^Y.h'' rfx' (x) (x') Vrs (x - x') * 

rs 

+ h.c.. (55) 

Here the kernel 



K.s(x) 



(27r)3 



(56) 



does not depend on the integration volume V and, under 
reasonable smoothness assumptions for Yrs,k, is a quite 
localized function around x = 0. 




empty space surrounding the atomic system. These are 
created inside the system as polaritons, and then trans- 
form into free-space photons when traversing the inter- 
face to the external vacuum. For the sake of simplicity, 
we shall assume that interface reflections are negligible 
at the frequencies of interest, so that the internal polari- 
ton is adiabatically transformed into the photon state at 
exactly the same energy. As light is collected far outside 
the slab, we do not have to consider here the case of light 
being guided inside the slab by total internal reflection. 

As usual in scattering problems, the final eigenstates are 
labelled by the wave vector k outside the atomic system 
and their frequency is fixed by the free-space dispersion 
Lo = ck. Outside the Mott insulator, polaritons reduce in 
fact to photons. Their wave function is the plane wave: 



V'r*(x) 



1 



(57) 



Inside the system, the same eigenstate corresponds to a 
polariton plane wave of the form [53| 



(58) 



The wavevector k™* and the branch index f are fixed by 
energy conservation in the refraction process at the inter- 
face separating the atomic system and the surrounding 
vacuum: the translational symmetry of the system along 
the X, y in-plane directions guarantees that the x, y com- 
ponents of the wave vector are conserved k^.y = fc™*, 
while the band index f and the z component fc™* are 
fixed by the energy conservation uifij^int = ck [58| . 

The amplitude A also depends on the external wavevec- 
tor k and is fixed by the particle flux conservation condi- 
tion. As interface reflections are assumed to be negligible, 
this reads: 



c cos(6') = 1^1 



2 Qr 



(59) 



FIG. 4: Sketch of the slab geometry under consideration. 



The advantage of the real-space Hamiltonian ([55]) is that 
it is not limited to bulk systems, but can be also applied 
to finite geometries. For the sake of simplicity, we shall 
limit ourselves to the plane-parallel slab geometry of figH) 
The size of the integration box in the z direction or- 
thogonal to the slab plane is assumed to be much longer 
than the slab thickness L^iab- In the transverse x,y di- 
rections, the slab is assumed to fill the whole integration 
box of size L,j: y. 

In order to apply the Fermi golden rule, one has to iden- 
tify the final states of the process. In the present case, 
they consist of pairs of radiative photons emitted in the 



where 9 (B™*) are the angles between the z axis and the 
propagation direction outside (inside) the slab: note the 
dramatic increase of the polariton amplitude inside the 
slab for v^ ^ c, a well-known effect in the theory of slow 
light propagation [3l| . Provided the slab is much thinner 
than the integration box ^ Lsiab, it is important to 
remind that the normalization and the density of states 
do not depend on the slab size, and only involve the total 
volume V = L^LyLz- 

Inserting the explicit form of the polariton wave func- 
tion ■0k(x, r) into the parametric Hamiltonian (|55p and 
using the Fermi golden rule, one immediately gets to an 
expression for the number of photons emitted into the ex- 
ternal vacuum per unit time dt, unit surface dT, = dx dy, 
and unit phase-space volume of transverse momentum 
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dN 



dtdT. d'^kj_ 

1 



(27r) 



L 



slab 



COS 6*^^^ (v^''^,. 



(60) 



For each value of the modulation frequency lo, the in- 
ternal wave vector A:*"* and the branch indices r, s are 
selected by the resonance condition (|^^ with k replaced 
by k™* . Inside the system the wave vectors of the emitted 
polariton pair differ in fact from perfect anti-parallelism 
by a negligible amount Ak^ oc l/Lgiat, still their fre- 
quencies W(j; j) j,i,it and their group velocities 
can be significantly different as soon as distinct branches 
are considered f ^ s. 

It is interesting to note that the result (I60|l is in perfect 
agreement with the emission rate predicted in (|T7)l for a 
bulk system of volume Vsiab = L^LyLsiab- in the absence 
of losses, all the polaritons created in the finite slab are 
in fact emitted from the system as radiation. 

By means of a change of variables, the emission rate ([HO]) 
can be rewritten in its final form as an emission rate per 
unit surface and unit solid angle dfl: 



dN. 



± 



1 



dt dfl 



(27r)2 
\ c 



V-, 



tint 



COSU± 



lab 



(61) 



where the ± index respectively refer to the photon which 
is emitted from the slab in the positive (-I-) or negative 
(— ) z direction at angles respectively 9± with the normal. 
This means e.g. that cj+ = ^i„t and uj^ — ojg -^int . 
The difference in the angular emission density in the ± 
directions is due to refraction effects of the f, s polaritons 
at the system- vacuum interfaces. 



A. Resonant enhancement in a cavity 

A possible way to further enhance the emission intensity 
is to surround the slab with a pair of mirrors of good 
reflectivity i? < 1: if the modulation frequency uj is on 
resonance with a pair of cavity modes, the dynamical 
Casimir emission results increased by a factor propor- 
tional to the finesse of the cavity [2l| . 

For simplicity, let us consider a plane-parallel mirror ge- 
ometry, so that the cavity modes are labeled by the in- 
plane component k^ of the wave vector, the branch index 
f, as well as by a positive integer number M defining the 
mode order along z. Assuming for simplicity that the 
mirrors are metallic and that no vacuum space is left 
between the slab and the mirrors, the polaritonic cavity 
modes are defined by the round-trip quantization condi- 
tion: 



slab 



ttM, 



(62) 



and their frequency is equal to the bulk polariton disper- 
sion fcint at the relevant wave vector /c*"* . Because of 
the non-trivial shape of the dispersion law ujf j^int as a 
function of fc™* , the cavity modes are not equally spaced 
in frequency. Their polaritonic wave function has a sim- 
ple sinusoidal form (the mirrors are at z = 0, Lgiab)'- 



ipf,MM^ (x, r) 



LslabLxLy \ Lslai 



(63) 

and their radiative decay rate into externally propagating 
photons is equal to 



1 - R 



1 = 



cos 



(6* )v^ 



gr 



(64) 



As expected, the slower the group velocity v^*", the 
smaller the decay rate 7. 

Fermi golden rule can again be used to estimate the emis- 
sion rate in the cavity geometry. As final states, pairs of 
cavity photons have to be considered, with a finite line 
width equal to 7 [s^. For a thick enough cavity, the 
mode order at the frequency of interest is M 3> 1 and 
overlap factor strongly privileges polariton emission into 
pairs of cavity modes of the same order M; the efficiency 
of all other processes is suppressed by their spatial phase 
mismatch. The emission rate into a pair of such modes 
is then easily obtained: 



dN 



dtdJ:d^k± 7r2 



It 



{cu~u;Tr + :f 



TT^ 7T 



(65) 



where ujt and jt are here the sum of respectively the 
frequencies and the linewidths of the pair of modes under 
examination. 

By comparing equation (|65p with the result (|60p in the 
absence of the enclosing cavity, and using ([M]) , it is im- 
mediate to see that the emission rate for an excitation 
exactly on resonance with a pair of modes {lu = lot) is 
enhanced by a factor 4/(1 — i?). This enhancement effect 
is quite general and holds for a variety of optical pro- 
cesses [11; for well reflecting mirrors i? < 1, it can be 
quite dramatic. 



VI. QUANTITATIVE DISCUSSION AND 
EXPERIMENTAL CONSIDERATIONS 



In the previous sections, we have obtained simple ana- 
lytical expressions relating the emission intensity in the 
different geometries to the Ks,fe parameter (|^ which 
carries informations on the microscopic optical proper- 
ties of the atomic medium. In the present section, we 
shall conclude the study by providing quantitative esti- 
mation of the emission intensity for realistic systems, and 
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discussing the most relevant issues that are hkely to arise 
in the design of an actual experiment. 




FIG. 5: Spectrum of the dynamical Casimir emission rate per 
unit surface and unit solid angle around the normal direction 
(^9± = 0) as a function of the modulation frequency uj for 
the resonant case uje = (jJm, a slab thickness Lsiab = 10 /im 
in the absence of enclosing cavity, and a relative modulation 
amplitude S^lc/^c = 0.05. Same system parameters as in 
FigEl (a) weak dressing Qc /C = 1.86-10"* < 1; (b,c) strong 
dressing flc/C — 2. The circles in (a,b,c) are the result of 
the exact calculation (|61|) . The solid lines in (a,b) are the 
analytical approximation (|50p . The red dashed line in (c) 
is the prediction of the dispersionless, time-varying dielectric 
model based on (|53|l . Panels (a,b) correspond to the band 
dispersion and absorption spectra shown in Fig[3ja-d). 



An example of emission intensity spectrum is shown in 
Figl5]for a Lgiab = 10 /im slab of Rb atomic Mott insula- 
tor in the absence of cavity as described by (pT|) . The cou- 
pling amplitude flc is assumed to have the periodic time- 
dependence (j49p at a frequency to, while its frequency loc 
is kept fixed: the modulation of then consists of a pair 
of coherent sidebands at ±2(a;e ~ ^g) around the carrier 
frequency uc — \(^e — w,„|. As expected, the central MP 
polariton branch appears as the most favorable region 
thanks to the combination of reduced group velocity and 
significant resonant mixing of light and matter excita- 
tions. In the resonant region [panels (a,b)] the result of 
the complete calculation is compared with the analyti- 
cal approximation (|50p : the agreement is excellent. In 
panel (c), the same calculation is performed in the low 
frequency region and is compared to Law's result ([55]) 
for a dispersionless medium [i^: inserting in (|53p the 
refractive index variation that follows from time depen- 
dence of r^c and using the general formula (pT|) for the 
emission rate, one obtains the dashed curve in Fig[S}:. As 



expected, at very low frequencies the agreement is excel- 
lent, but dispersion effects start playing a significant role 
already for lu « 0.2 cje. 




10"" 10"' lo" 10' 10^ lo' 



n„/co 

C e 

FIG. 6: Prediction (|6ip for the dynamical Casimir emission 
rate per unit surface and unit solid angle around the normal 
direction {6± = 0) as a function of the coupling amplitude 
Sic for a fixed value of the relative modulation amplitude 
5Qc/^c ~ 0.05 at the fully resonant point uj/2 = lue — uj,n. 
Same system parameters as in Fig[5] the magenta circle and 
the blue diamond correspond to respectively the weak (a) and 
strong (b) dressing regimes. 



All panels of Fig[5] have been calculated for the same 
value of the modulation amplitude SQc/^c- in spite of 
the higher value of v^'', a larger value of fie is favourable 
in view of maximizing the emission intensity as it allows 
for a stronger modulation of the optical properties. A 
specific plot of the resonance emission rate as a function 
of the coupling amplitude flc is shown in FiglGl as pre- 
dicted in Sec lIV Bl the Casimir emission first grows as il^ 
for flc <C C, then saturates to a finite value for flc 3> C. 

Remarkably, the emission intensity at resonance can 
reach quite substantial values already in the absence of 
a cavity. In the weak dressing case {ftc/'^Tr — 12 MHz), 
the rate of emitted photons from a 1 cm^ system in the 
unit solid angle around the normal is of the order of 1 
photon per second [see FiglHJa)]. Then it quadratically 
(see Figini) increases for growing V,c, to eventually sat- 
urate around a value larger than 10^ photon per second 
for a huge dressing ric'/27r ~ 100 GHz [see FiglSjb)]. For 
alkali atoms such as Rb, a dressing amplitude flc/'^n in 
the 10 MHz range corresponds to ijntensities of the dress- 
ing beam in the mW/cm^ range |i39j . 

A crucial difficulty of most dynamical Casimir experi- 
ments consists of varying the optical properties of the 
system at a high enough speed. In our specific setup, 
this amounts to modulating the dressing beam at a fre- 
quency which is resonant with the creation of a pair of 
MP. 

Although this is hardly done with the almost symmetric 
A schemes currently used in slow light experiments with 
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alkali atom samples [T^, [H, [3| , still strongly asymmetric 
A schemes as shown in FigHJc) can be used. The ground 
state of the atom is now the m state, while the g state is 
a long-lived, high-energy metastable state. The dressing 
beam then acts on the to — > e transition at a frequency 
jo^e — Wm| higher than twice the frequency jcje — cjgj of the 
(7 — > e transition on which the dynamical Casimir radia- 
tion is to be emitted: no principle difficulties then appear 
to prevent one from modulating the dressing amplitude 
ric at the required frequency w ~ 2 jwe — w^j [59|. 
Suitable strongly asymmetric A configurations can be 
found e.g. in alkali-earth atoms [i^ whose laser cool- 
ing and trapping techniques have experienced remarkable 
advances in the last few years A specific choice in 
view of our dynamical Casimir application can be ^^Sr 
atoms, whose 4d^D2 metastable state appears to have 
the required properties to be used as the g state of a 
strongly asymmetric A scheme. It is in fact connected to 
the excited bp^Pi state by an optically active infrared 
transition at A = 6.5 fim which, for a lattice spacing of 
300nm, gives a value C/oJe — 5 • 10~^ not far from the 
Rb one. The Sp^Fi state can be dressed by driving the 
atom on the A = 461 nm transition from the absolute 
ground state of the atom 5s^ -^Sq 5p^Pi. Before per- 
forming the dynamical Casimir experiment, atoms have 
to be optically pumped in the Ad^D2 metastable state 
e.g. by means of a tt Raman pulse; the lifetime of the 
state being of 0.33ms, there is enough time left to carry 
out the dynamical Casimir experiment before atoms de- 
cay to the ground state via the 5p^Pj state: even for a 
group velocity as low as lOm/s, the transit time across 
a 10 ^m thick cloud is in fact much shorter, of the order 
of 1 ^s. Note also that this spontaneous decay channel 
involves photons at 1.8 /im and 689 nm, in completely dif- 
ferent spectral regions from the dynamical Casimir ones 
at 6.5 /im, which can therefore be spectrally isolated with 
no difficulty. 

Another, may-be simpler solution is to stick to alkali 
atoms such as Rb, but use a ladder scheme [see Fig[T][b)] 
where the to state is a electronically highly excited state 
instead of a A one. Even though the unavoidable ra- 
diative contribution to "f„i prevents the EIT effect from 
completely killing the absorption, still DCE light can es- 
cape the slab without dramatic losses. This, at least for 
the large values of flc of present interest [FigHI^d)]. 

As the EIT effect is not crucial for the observation of 
the dynamical Casimir emission, one may think of do- 
ing experiment with atomic samples where spontaneous 



emission from both the e and the to states is still ef- 
fective, e.g. Bose-Einstein condensates or even thermal 
gases. Additional difficulties could however arise in such 
systems from the inhomogeneous broadening of ujm due 
to the spatial variations of the trapping and interaction 
potentials and from the reduced value of the light-matter 
coupling coefficient C due to the reduced atomic density. 
While this latter effect can be quite dramatic in non- 
degenerate clouds, the inhomogeneous broadening of Um 
should be in both cases overcome by choosing a suffi- 
ciently strong coupling amplitude flc- 



VII. CONCLUSIONS 

In this paper we have performed a complete systematic 
analysis of the optical properties of a gas of coherently 
dressed three-level atoms trapped in an optical lattice in 
a Mott insulator state. The extreme degree of coherence 
of this system allows for propagation of light at ultra-slow 
group velocities for long times and distances. The optical 
properties of the medium can be controlled in real time 
by varying the amplitude and the frequency of the dress- 
ing field. For sufficiently fast modulation rates, the zero- 
point fluctuations of the polariton vacuum state are con- 
verted into observable radiation by dynamical Casimir 
effect. 

We have developed a general theory to quantitatively 
characterize the dynamical Casimir emission in terms of 
a simple parametric Hamiltonian and we have identified 
the most favourable case of a resonant dressing frequency 
whose amplitude is periodically modulated in time. Ex- 
perimentally realistic geometries such as plane-parallel 
slabs and planar cavities are analysed in detail. Remark- 
ably, a sizeable radiation intensity is predicted for state- 
of-the-art systems and no spurious emission from black- 
body radiation or incoherent luminescence is expected to 
mask the dynamical Casimir signal. 
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The commutator of raising = \^)j{g\j Bind lower- 
ing ae,j = \g)j{e\j operators is [ae,j,alj = \g)j{g\j - 
= Ij — \''n)j{m\j — 2 \e)j{e\j, where Ij is the iden- 
tity operator on the jth atom. If the occupation of both 
e and m states is much smaller than 1, this commuta- 
tor is well approximated by the identity Ij and thus has 
bosonic nature. 

Coupling to the longitudinal e.m. field corresponds to 
the static dipole-dipole Coulomb interaction [^, H^]. At 
low excitation regimes, this term is only responsible for 
a slight red-shift of the g —> e transition frequency uie 
by AtJe = —2C^/{Suje). For the typical values of the 
systems under examination here, this shift is on the or- 
der of a few 10 MHz and can be reincorporated in the 
definition of oje. In a semi-classical description of light- 
matter interaction, such a shift naturally appears when 
the Clausius-Mossotti form of the dielectric constant is 
used [33, 

Though small, this term is crucial to preserve gauge in- 
variance and avoid unphysical behaviors in the fc — > 
limit. 

In our formalism, such processes correspond to the con- 
version of a e excitation into a m excitation plus a po- 
lariton at a frequency u>e — Wm far from the resonance 
region uje — (^g- Using Fermi golden rule, it is immediate 
to see that the corresponding decay rate is indeed equal 
to the free space spontaneous radiative decay rate for the 



e ^ m transition. 

[54] Note that this plot does not include the extinction effects 
that occur for frequencies within the (very small) gaps 
between the bands, and which are e.g. responsible for 
the operation of distributed Bragg reflectors [H. |4^. 

[55] For the Rb atoms under consideration here, an absorp- 
tion length of 1 m corresponds to k"'^" /ke — 1.24 ■ 10^^. 

[56] A similar perturbative calculation for a periodic modula- 
tion of uim would lead to a vanishing result at the lowest 
order in perturbation theory. This is a direct consequence 
of the vanishing value of u™^^ obtained in p4[) . A com- 
plete calculation including next order terms results in a 
dynamical Casimir emission intensity orders of magni- 
tude weaker than for a modulation of Q,c- 

[57] In real space, the polariton wave function i/)Jj"(x, r) de- 
pends in fact on the spatial coordinate x and on the 
branch index r. This latter can be interpreted as a sort 
of internal degree of freedom like spin. 

[58] As the bands have no spectral overlap, at most one prop- 
agating polariton mode exists for a given frequency oj: 
this guarantees that no additional boundary conditions 
are required here [13]. Rather, the energy conservation 
condition can be satisfied only outside the (very small) 
gaps that open in the polariton dispersion between the 
LP and the MP, and between the MP and the UP. For 
frequencies inside these gaps, the radiative eigenstates 
are confined in the external vacuum and can not pene- 
trate the atomic system: the dynamical Casimir emission 
is therefore strongly suppressed. 

On the other hand, light escape from the atomic system is 
possible only for polariton modes which lie within the so- 
called radiative cone + < uP /c?: polaritons which 
are created outside this cone remain in fact trapped in 
the slab and propagate through it as in a waveguide. 
[59] A dressing amplitude Q.c{t) modulated at the optical fre- 
quency uj can be obtained e.g. by nonlinear wave mixing 
on a suitable nonlinear crystal of a carrier at ujc with 
another beam at uj. 



